Surface Integral of Babu Diagram 
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We point out that the surface integral of the Babu diagram gives non-trivial main contribution 
independent of the charged lepton masses and the scale of new physics. It can fit the experimental 
data on both solar and atmospheric neutrino oscillations. However, then the coupling constants of 
the new Higgs singlets with the leptons gain values smaller than that in the usual two-loop radiative 
analyses. This yields distinction between our consideration on ultraviolet-behaviors and the former 
based on the infrared radiative corrections. 

PACS numbers: 12.60.-i, 12.15.Lk, 14.60.Pq 
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The evidence of neutrino oscillation experimentally in- 
dicates that neutrinos are massive particles and that fla- 
vor lepton number is not conserved. Since in the stan- 
dard model (SM), neutrinos are massless and flavor lep- 
ton number is conserved. The neutrino oscillation exper- 
iments are a clear sign that the SM has to be extended. 

There existed two different typical approaching ways 
on generating neutrino masses such as the seesaw mech- 
anism (SSM) Q and the radiative mechanism (RM) 
0. The SSM has been the most popular explanation 
of the small neutrino masses. It is simple, and re- 
lies only on dimensional analysis for the required new 
physics. Current neutrino data point to a seesaw scale of 
Mr ~ 10 10 — 10 15 GeV, where lepton number violation 
occurs through the Majorana masses of the right-handed 
neutrinos. With such a high scale, the effects of lepton 
flavor violation in processes other than neutrino oscil- 
lation itself become extremely small. For example, the 
branching ratio for the decay n — » e + 7 is of order 10~ 50 
within the seesaw extension of the SM. 

An alternative to the SSM, the RM also explains the 
smallness of neutrino masses naturally. In this approach, 
neutrino masses are zero at the tree level and are induced 
only as finite radiative corrections. These radiative cor- 
rections are typically proportional to the square of the 
charged lepton (or quark) masses divided by the scale 
of new physics. In this case, lepton flavor violation in 
processes other than neutrino oscillations may become 
experimentally accessible. 

On RM, Babu has introduced a model Q (called Babu 
model) which is directly extended from the SM by thank- 
ing to two charged singlet Higgs fields to induce neutrino 
masses at the two-loop level. It has been judged that 
the neutrino masses arise to be naturally small. How- 
ever, if one takes to carefully calculating that two-loop 
diagram (called Babu diagram) , the one will see it yields 
additionally a radiative mass term which is much larger 
than the term derived by his comments and calculations. 
Moreover, this term completely results from the ultravi- 



olet behaviors with large internal momenta of the dia- 
gram. Further, if taking evaluation for neutrino masses 
and mixing, it can fit not only experiment data but also 
yields difference from that based on the usual radiative 
corrections at low energy. 

The aim of this Letter is to clear up these important 
remarks and show some interesting consequences. 

The Babu model Q includes two SU(2)i singlet Higgs 
fields, a singly charged field h + and a doubly charged 
field k ++ . Moreover, right-handed neutrinos are not in- 
troduced. The addition of these singlets gives rise to the 
Yukawa couplings: 



Cy = fab(lpaL) C 1pbLh + + h ab (l aR ) C l bR k ++ + Jl.C, (1) 

where ipL stands for the left-handed lepton doublet, Ir 
for the right-handed charged lepton singlet and (a, b) 
being the generation indices, a superscript c indicates 

charge conjugation. Here ip c = Cip with C is the 
charge-conjugation matrix. The coupling constant f a b 
is antisymmetric (f a b = —fba), whereas h a b is symmet- 
ric ih a b = hba). Gauge invariance precludes the singlet 
Higgs fields from coupling to the quarks. Eq. (JIJ con- 
serves lepton number, therefore, itself cannot be respon- 
sible for neutrino mass generation. The Higgs potential 
contains the terms: 



V(<f>, h+,k ++ ) = [i{h,-h,-k ++ + h+h+k—) + 



(2) 



which violate lepton number by two units. They are ex- 
pected to cause Majorana neutrino masses. 

On Babu's opinion in his model, Majorana neutrino 
masses are generated at the two-loop level via the dia- 
gram shown in [j| and again depicted in Fig^ Then, he 
has written down the corresponding mass matrix element 
for Majorana neutrinos as follows 



Snf a ch cd m c mdlcd(f + )db, 



(3) 



in which h a b = rjh a b with 77 = 1 for a = b and n = 2 for 
a ^= b. The integral I cc i is given by 

d A k f d 4 q 
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FIG. 1: The two- loop diagram in the Babu model. 



mental data 0, Q . 

However, now we will point out that the above solution 
of M a i, is not complete. To see this, we will carefully 
apply the Feynman rules to the diagram. After that, the 
large corrective term is followed which kept in his model. 

In the momentum space, the propagator of the charge- 
conjugate spinor field is given as follows [j| 

(0|T(^ c V^)|0) = C{Q\T{^)\Q) T C- 1 
-^CS T (p)C- 1 = {-j>- m y l =S(-p). 



After evaluating |0J, it has been seen that the neutrino 
masses are small finite and compatible with the experi- 



For simplicity, f a b and h a b are assumed to be real. Then, 
the mass matrix of the Majorana neutrinos which is in- 
duced by the above diagram reads 



iM' 



'^(-S^fachcdfdb) / 



d 4 k d 4 q 



P L (Jz+ m c )P L {-fr + m d )P L 



(2tt) 4 (2tt) 4 (fc 2 - m 2 c ){q 2 - m 2 d )[{q - k) 2 - m 2 ](fc 2 - m 2 ){q 2 - to 2 ) ' 



(5) 



The P^s in the vertex factors have to be kept as coeffi- 
cients out of the partial momentum integrals defined for 
the free spinor and boson fields. Even if their overall inte- 
gral is obtained, they still are in the coefficient. Keep our 
mind that, here the © has conveniently been written. 
We have 



P L $+m c )P L (- 



m d )PL = 

Pl^PJPl + PL^ L m d 
m c P L iPL + m c m d P L . (6) 



Therefore, the r.h.s of J^J is summed of four terms re- 
spectively to the terms in the r.h.s of ||BJ|. The fourth 
is finite which has been considered by Babu in his pa- 
per. Thus, applying its Pl on the left-handed Majorana 



field in the mass Lagrangian, e.g. — ^{vl) c M'^ A > P^v^ — 
-\{v L ) c M'^v L , we get then: 



M 



'(4) 



= 8^ 2_j f ached fdbm c m d I cdl 



(7) 



with I cd given in For the second and the third, their 
integrals are finite, while the respective coefficients are 
vanish. Thus they are explicitly canceled: 



M. 



'(2) 



M. 



'(3) 



0. 



(8) 



Finally, to keep the PlS in the coefficient of the re- 
spective overall momentum integral, the first term can 
be explicitly given in the following form: 



-iM'^ = 8iflJ2fachcdfdbPL7t,PL-/uPL f 



d 4 k d 4 q 



k^q" 



(2tt) 4 (2tt) 4 (fc 2 - m 2 c ){q 2 - m 2 )[{q - k) 2 - m 2 ](fc 2 - m 2 )(q 2 - m 2 ) ' 



(9) 



It can be checked that the two-loop integral of fc, q in the 
r.h.s of © is logarithmically divergent while its coeffi- 
cient is vanish, hence 



- iM ab ] = CXOO, 



(10) 



which follows an indeterminate radiative term of the 
Babu diagram. 

Explicitly, this term has been canceled in the former 



analyses Q which present a fallacy. Suppose 



PLhPUPL = P L Pl 



(11) 



happens in JSJ. Indeed, the source of the fallacy is that 
algebraically permuting of the matrices such as Hll|) in 
the integrand of the logarithmically divergent momentum 
integral of (JSJl changes the value of the integral by an 
finite amount, which is equal to l|10|) (also see 0). 
To evaluate (|10|) , dimensional regularization can be de- 



3 



veloped. Since the anticommutativity of 75 in d dimen- 
sions is lost Q (also see ||), the idea is that [To||: The 
external indices and the momenta all live in the phys- 
ical four dimensions, the loop momenta k, q have com- 
ponents in all dimensions. 75 is expanded continuously 
analytic to d-dimensions such that anti-commutes with 
7^ for fj, = 0, 1, 2, 3 but commutes with 7^ for other val- 
ues of ll. Let k = fc|| + k±, q = q\\ + q± we have then 

m 



Pl^Pl = P L h 



Pt, 



-kq. 



(12) 



The coefficient in <|12|) owns d — 4 which yields zero-limit 
whenever returned back to the four dimensions, hence we 
need only to define the divergent parts of the integral in 
@. By denoting 4 — d = 2e, it can be checked that 



/ = 



d d k d d q 



kq 



(2ir) d (2ir) d (k 2 
1 



m c)(9 2 — rn d)[( < l ~ 
1 



(k 2 - m 2 )(q 2 ~ m 2 ) 2(47r) 4 e 
Then, the indeterminate term yields 

x 00 = 8ifi y~] f ached fdbPh —j— x I 



0(1) 



(13) 



c,d 

-Aifl ^ / gjlcd fdbPh I < 
c,rf 

-^^V^ f ached fdbPL + 0(e) 



2(4tt) 



0(1) 
(14) 



By shifting the Pl into the left-hand Majorana field in 
its mass Lagrangian, we can rewrite 



M 



'(i) _ 



(47T) 



X)/«cW«tt + 0(e)- (15) 



Hence, with the help of Q, © an d CSJ, the mass matrix 
of the Majorana neutrinos in the four dimensions is given 
as follows 



(4tt)' 



[M /] fachedfd 



8/i fach c dfdb i mcm d I cd . 
c . d 



(16) 



The following remarks are in order: 

1. The radiative term (jl5(l is dependent only on the 
Yukawa coupling constants of new physics such as 
fab, h a b, but not on the masses of the charged lep- 
tons and of the new Higgs scalars. Thus, it should 
also be present in the massless theory. 

2. It can be checked that this term is not affected by 
higher-order radiative corrections, e.g. the Babu 
diagram with more than two loops do not con- 
tribute to this term. 



3. The term completely results from the ultraviolet 
asymptotic (large momenta) of the integrand in © 
and thus in ©, is also identified as the surface 
integral of the diagram. Let us analyze that, the 
integrand in the Q is strongly local, large around 
zero-momenta and hence gives main distribution to 
this type of the RM, as usual (infrared behavior); 
while the (|15|l is distinct from that upon ultraviolet 
behavior. 

4. On the evaluation of the Babu integral I c d in Q 
(also see 0, 0, Ej), the terms in JjJ are propor- 
tional to, e.g. m c irid/A with A ~ m^m 2 ,, there- 
fore, it is very smaller than that in the (|15f) . For ex- 
ample, choosing ~ m/, ~ 100 GeV, the largest 
one is only m 2 /K w 0.00031 < 1. Thus the surface 
integral gives 771am contribution on the Babu dia- 
gram. The mass matrix <|16|) can be approximated 
by 



ML 



(An) 



'^y^ fachcdfdb, 



(17) 



which completely differs from early that in © . 

It can be checked that the matrix M' ab in (|16|l as well 
as its approximation in l|17|l always has a zero eigenvalue 
independent on its parameters; and the respective eigen- 
state is dependent only on the /„;, couplings, which is 
identified as: 



m Vl = 0, 



\j fefi fer fjlr 



(18) 

(f^rVe ~ fer^p, + fe^r)- 

(19) 



Next, from the orthogonal condition of the v\ eigenstate 
to the two remaining neutrino vectors, one can write 



1 



i er f lit 



-(ferVe + Ut^/j.)} 



--(fe^ + ferVr). (20) 



Then in this basis of (v\, v^v's), the matrix M' ab in lltj|) 
owns the form: 




]j" = | /»,..; m v > 2U ' 3 
m v > v i m v > 



(21) 



The matrix M" fc gives mixing of v' 2 and v' 3 . Next, by 
rotating over the mixing angle (f> in the plane (v' 2 , u' 3 ) — ► 
(i>2, z/ 3 ), we have 



tan 2(f) 



2m v 



(22) 
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and the physical mass eigenvalues are given by 



(EU yield 



m v i + m v i — \ (m u i — m v i ) 2 + 4m 2 , , 



m„> + m,.i + , / (m,/ — m,y ) 2 + 4m 2 , , 



(23) 



(24) 



With the help of i|19|) , H2U|) and the <fi, we can easily ob- 
tain the mixing angles of the three neutrino generations 
in this model respective to the standard parametrization 
of an 3 x 3 Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
[l2l IT3II unitary mixing matrix: 



in. 



2fi 



[(fl+f'r)h T7 



f 2 

J ey 1(2 u 1 (2 



f2 , f2 
J er 1 J lit 



+ 



(4tt; 
f 2 

J e^i ' J er 



Ml f f2 u , f 2 . \ 

Ue/i" TT ' Jer n ^n) 



feu, fur 



(31) 



(32) 



(My(/;+^)^+/, 2 T ) 

X [feri^uu — ^ ee — ^tt) — feu^ee ~ fur^rr]- 



(33) 



sinf 



tan( 



13 



/f2 1 f2 
\f J er J [it 



sin 0, 



(25) 



fer\/ fefj, + /er + //IT 



■ COS( 



//iT \ fer "T" fur 



(26) 



tan 6> 23 = 



//XT + /er Sin + /ep J fer + fjlr C0S < 



fer J fir + / 2 T COS ( 



(27) 



The equations (|28() and (|33|) give 



heeifL + fl) + Kr{fl + flr ) 



/, 



(34) 



In addition, since and are not mixed, the l|31|) and 
(El yield 



in,, 



2M (fL + fer + fl) Urr + h ^ 



(47T) 



f 2 + f 2 / ' 



(35) 



The equations H25|) . (|26[l and (|27l) yield that the mixing 
angles #12, $23 and #13 can be given by studying the <f> 
on the mixing between v' 2 and t^. Moreover, since the 
mixing angle #13 is small, e. g. t he CHOOZ experiment 
gives sin 2 #13 < 0.08 (also see H3), thus Eq. l(2*5|) requires 
that the mixing angle <j) is respective small too. 

Now we will show that above solution can fit the cur- 
rent data on both solar and atmospheric oscillations. For 
simplicity, we shall work with the matrix (|17|l by taking 
a special case of #13 = 0. It yields = 0, and then: 



= 0, 



tan #1 



tan B-) 



fer \ feu + fer + fur 



fury fer + far 



fen 
fer 



(28) 
(29) 

(30) 



Next, on the constraints in the various parameters 0. Q|. 
it has been agreed that ~ for a ^ d, then with the 
help of ljT7|) . (jT§|) and (|2"U|) the elements of the matrix 



1711/3 ~ (47r) 4 + ^ T + ^" T ^ ^y ee + hrT f 2 + f 2 ^ ' 

"(SO) 

Upon l|29|l and (|30|l . we see that #12 and 023 are not 

necessarily small. Then to fit with observations (also see 
13]), one can choose: 



tan 2 6*23 ~ 1, tan 2 6 12 ^ 0.43, 



which yield 



fep, — fer ~ 0.585/^ r . 



(37) 



(38) 



Denoting Am| = m 2 . — m 2 ^ , from experimental data 
[I3J on solar and atmospheric neutrino oscillations, 
one gets the following relation Am^/Am^ as 6.9 x 
10" 5 eV 2 /2.6 x 10~ 3 eV 2 . With the help of the (29, 
the equations l|33)l and (|3*^f) give the following approxi- 
mation: 



.123h TT , 



which yields 



1.782 X 10- 5 ^f 2 T hee, 

1.108 x lQ- 4 vf* T h ee . 



(39) 

(40) 
(41) 
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Let us put far — h ee |l| , then it follows 



— ^ee 



7.754 eV 



1/3 



/' 



1/3 



(42) 



Choosing /i = 200 GeV, which is the highest scale in 
the theory, we see that f/j, T ,h ee ks 1.33 x 10~ 3 . On the 
constraints (|54l . i|3~5|) and the remaining couplings 
gain: 



V » 2.01 x 10 



0.775 x 10~ 3 , 
3 h » - 



(43) 

1.63 x 10" 4 . (44) 



It is easy to check that all above values for the couplings 
explicitly satisfy the experimental constraints given in 
Moreover, the above solution can be also used to 
explain the mass hierarchy and mixing of Majorana neu- 
trinos at some energy scale wherein the gauge symmetry 
of the SM vacuum is not broken. 

Thus, we have derived a new main radiative term of the 
Babu diagram which can fit the data on both solar and 
atmospheric neutrino oscillations. However the required 
couplings /, h for the model are smaller than that in the 



usual two-loop radiative analyses (also see |l4||). The 
reason lies behind ultraviolet behavior on the radiative 
mechanism of generating neutrino masses. 

To conclude, in this Letter we have shown that 
the Babu diagram owns a non-trivial main term in 
ultraviolet-behavior. Its contribution is given as a nat- 
ural consequence of the surface integral which is defined 
in the term of dimensional regularization with the new 
properties of 75 in d-dimensions. After that, we have 
taken the exact diagonalization of the neutrino mass ma- 
trix with the help of the intermediate mixing angle <f>. 
Requirement of (f>'s smallness naturally gives the mixing 
parts consistent with current experimental data on both 
solar and atmospheric neutrino oscillations. This con- 
sideration predicts the coupling constants f a b and h aa 
are smaller than that in the former two-loop analyses 

It is also evident that the ultraviolet-asymptotic be- 
havior in the radiative mass mechanism should be men- 
tioned in the current neutrino problems. This conclusion 
adds one more nice feature to the Babu model. 

The authors wish to thank T. Kitabayashi and M. 
Yasue for reading the manuscript and comments 
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